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Small modeling errors in a finite-element model will eventually induce errors in the structural flexibility and mass, 
thus propagating into unpredictable errors in the unsteady aerodynamics and the control law design. One of the 
primary objectives of the X-56A Multi-Utility Technology Testbed aircraft is the flight demonstration of active flutter 
suppression and, therefore, in this study, the identification of the primary and secondary modes for the structural 
model tuning based on the flutter analysis of the X-56A aircraft. The ground-vibration test-validated structural 
dynamic finite-element model of the X-56A aircraft is created in this study. The structural dynamic finite-element 
model of the X-56A aircraft is improved using a model-tuning tool. In this study, two different weight configurations 
of the X-56A aircraft have been improved in a single optimization run. Frequency and the cross-orthogonality 
(mode shape) matrix were the primary focus for improvement, whereas other properties such as c.g. location, total 
weight, and off-diagonal terms of the mass orthogonality matrix were used as constraints. The end result was an 
improved structural dynamic finite-element model configuration for the X-56A aircraft. Improved frequencies and 
mode shapes in this study increased average flutter speeds of the X-56A aircraft by 7.6% compared to the baseline 
model. 


Nomenclature 

A = area of bar cross section 

E = Young’s modulus 

/ = frequency 

G = shear modulus 

1 1 - area moment of inertia of bar with respect to plane 2 

1 2 = area moment of inertia of bar with respect to plane 1 

J = objective function 

J k = performance indices 

K 2 = spring constant in direction 2 

K 3 = spring constant in direction 3 

K 4 = spring constant in direction 4 

M a = analytical mass matrix obtained from Nastran 

M = orthogonal matrix 

MAC = modal assurance criteria 

m = number of modes to be use in system equivalent 

reduction expansion process 
n - number of modes for model tuning 

n AD = number of A -set degrees of freedom 

n MD = number of measured (or master) degrees of freedom 

S = cross-orthogonality matrix 

V = velocity 

0 Aset - eigenmatrix obtained from Nastran modal analysis 

= intermediate matrix to compute M A 
O m = eigenmatrix, corresponds to master degrees of freedom 

O s = eigenmatrix, corresponds to slave degrees of freedom 


I. Introduction 

O NE of the major goals of the Fundamental Aeronautics program 
under the NASA Aeronautics Research Mission Directorate 
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(ARMD) is to develop a cutting-edge technology for higher- 
performance lighter-weight aircraft. Higher performance includes 
energy efficiency and operability technologies that enable advanced 
airframe and engine systems. Removing weight from an aircraft 
usually results in reduced stiffness and, therefore, increased flex- 
ibility. The increased flexibility creates an aircraft that is more 
susceptible to aeroelastic phenomena such as flutter, divergence, 
buzz, buffet, and gust response. Therefore, structural weight re- 
duction is constrained by these adverse aeroelastic phenomena. 

The ability for flexible motion control, such as gust-load 
alleviation and active flutter suppression, will help in opening the 
new design space for lighter- weight aircraft, which will increase fuel 
efficiency and help with noise reduction. Many active aeroelastic 
controller designs have been developed and proven analytically. In 
some cases, wind-tunnel tests have been performed to demonstrate 
these designed active controllers. Based on the maturation of active 
aeroelastic controls, it is time to perform an actual flight dem- 
onstration of these active controllers. However, there has been limited 
testing on active flutter suppression technology due to the high risk of 
vehicle loss. 

An open, nonpropriety, unmanned aerial vehicle called the Multi- 
Utility Technology Testbed (X-56A), as shown in Fig. 1, has been 
developed by the U.S. Air Force Research Laboratory (AFRL) for 
relatively high-risk aeroelastic flight demonstrations and future 
research. If a flight goes beyond flutter boundaries and a wing fails, 
then the X-56A aircraft will be fitted with a center-body-mounted 
ballistic parachute recovery system. Sharing the same goal to design a 
lighter-weight aircraft by suppressing associated aeroelastic insta- 
bilities using an active control has led to a joint collaboration between 
NASA and the AFRL. The AFRL provided the X-56A aircraft, 
ground control stations, analytical models, and ground and flight-test 
data to NASA. After the first flight test with flexible wings, AFRL’s 
assets, which include two reusable center bodies, one rigid wing, 
three flexible wings, and ground control stations, will stay at NASA 
Armstrong Flight Research Center (AFRC) for future technology 
demonstrations, such as the NASA-designed lighter-weight wing 
with curvilinear spars and ribs, the Northrop Grumman-designed 
flying wing, and/or the Boeing-designed joined wing, as shown in 
Fig. 2. Furthermore, Lockheed Martin Skunk Works (LMSW) has 
provided NASA the finite-element models and the ground-vibration 
test data for this research project. 

The X-56A aircraft has a wingspan of 28 ft and a weight between 
200 and 480 lb, depending on its configuration [!_]. This aircraft is 
powered by twin JetCat P400 turbojet engines and configured for 
easy wing replacement. The design also includes a hard point on the 
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Fig. 1 X-56A aircraft with ground control stations. (Figure courtesy AFRL/Lockheed Martin). 



Fig. 2 Future research configurations of the X-56A aircraft. 


center upper deck of the aft fuselage that can be adapted to house 
either a third engine or the boom for a joined wing, thereby enabling 
testing of more advanced aerodynamic concepts. A ground- vibration 
test was completed on 17 December 2012, and test data have been 
provided to NASA for model validations. 

Modem aircraft make extensive use of composite materials to 
reduce weight. Aircraft aeroservoelastic models have significant 
levels of model parameter uncertainty due to the composite manu- 
facturing process. In an attempt to mitigate this, ground- vibration test 
(GVT) data can be used to tune a model and help to improve its 
accuracy over a wider range of operating environments. This tuning 
and adaption procedure has become an accepted and necessary step in 
making models more reliable; however, it cannot handle or eliminate 
all uncertainties in a finite-element (FE) model. To minimize the 
uncertainties in aeroelastic as well as aeroservoelastic response com- 
putations, the stmctural dynamic model, the unsteady aerodynamic 
model, and the actuator model should be validated with respect to the 
corresponding test data [2,3]. Small modeling errors in an FE model 
will eventually induce errors in the structural flexibility and mass, 
thus propagating into unpredictable errors in the unsteady aero- 
dynamics and the control law design. Tuning an FE model using 
measured data to minimize the modeling uncertainties is a chal- 
lenging task in the area of structural dynamics. 

Over the past few decades, major progress has been made to 
develop methods to update these finite-element models to match 
computed frequencies and mode shapes to the corresponding 
measured GVT data [4-8]. Recently, attention has shifted toward 
applying these methods to the health monitoring of structures to 
address this problem [9-12]. Literature reviews on finite-element 


model updating methods are detailed in [13-15]. These methods can 
be separated into two different categories: direct methods and indirect 
(or parametric) methods. The direct methods correct the mass and 
stiffness matrices but do not take into account the structure’s physical 
characteristics. On the other hand, the indirect (or parametric) 
methods correct the model by altering the physical structural pa- 
rameters, which include (but are not limited to) structural sizing 
properties, point mass and spring properties, and material properties. 
In this study, the parametric method was selected for the model- 
tuning optimization process. 

Supporting the ARMD guidelines, NASA AFRC has developed a 
computer code for FE model tuning [2,16] using the object-oriented 
optimization (O 3 ) tool [ 17 ] together with Nastran [18], an FE analysis 
computer software program. Three optimizer algorithms have been 
incorporated into the O 3 tool: design optimization tools (DOTs) [19], 
based on a gradient-based algorithm; the genetic algorithm (GA) 
[20], which is a class of stochastic and global optimization; and big- 
bang-big-crunch algorithms [21-24]. This O 3 tool, as shown in 
Fig. 3, was crucial in the development of a multidisciplinary design, 
analysis, and optimization tool [ 25 ] and an unsteady aerodynamic 
model-tuning tool [3]. 

The primary objective of this study is to obtain the GVT- validated 
structural dynamic finite-element model of the X-56A aircraft. Based 
on the validated model, flutter boundaries are computed again and 
compared with results obtained from the nonvalidated model. The 
secondary objective of this study is to improve and extend the 
previously developed structural dynamic FE model-tuning tool [2] 
and demonstrate this tool using the X-56A aircraft as a sample 
case. 
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Fig. 3 Structural dynamic model tuning based on the O 3 tool. 


II. Optimization Theoretical Background 

In the FE model-tuning tool, the user chooses an optimization 
algorithm, defines the continuous as well as the discrete design 
variables with side constraints, starting values, and scaling factors, 
and provides script and data file names for each performance index. 
Communication between the O 3 tool and each analysis submodule is 
performed through the use of design variables and performance 
indices. Once the O 3 tool is executed, design variables are created and 
saved in an external data file. This data file cannot be shared with 
other analysis submodules; therefore, a script command makes a 
copy of this data file first. Script commands for each performance 
index perform the following tasks. 

1) Prepare data for analysis program (preprocessor module). 

2) Submit analysis job (analyzer module). 

3) Read result files, compute the required performance index 
value, and save it (postprocessor module). 

Analysis submodules of the extended model-tuning tool are 
explained in this section. Blocks in Fig. 3 can be categorized into 
three modules: the preprocessor, the analyzer, and the postprocessor 
modules. The preprocessor module (update Nastran input deck 
block) in Fig. 3 is used to create and update input files based on the 
design variable values provided by the O 3 tool before executing the 
analyzer module. The analyzer module (Nastran modal analysis 
block) in Fig. 3 can be a commercial and/or the in-house codes for a 
specific discipline. The script commands will execute the analyzer 
module automatically. The postprocessor module is used to post- 
process the Nastran output file, which is computed from the analyzer 
module, and it automatically computes the performance indices. 
Several submodules were developed for the structural dynamic FE 
model tuning. 

A. Weight Submodule 

A weight generator in Nastran code and measured test data are used 
to compute the maximum 10 performance indices associated with the 
total weight, c.g. locations, and the moments of inertia. Note that the 
G in the subscripts for the following variables displayed represents 
the test-measured data value as shown in Eqs. (1-10): 

(W — W r \ 2 

h = ( ^ ) (total weight) (1) 

( X — X \ 2 

x - ) (x c.g. location) (2) 


H 

Y — Y g \ 2 

Y J O' c-g- location) 

(3) 

h = ( 

'Z - Z G \ 2 

— — 1 (zc.g. location) 

(4) 

J S = 

(I XX lx xG 

V hxG J 

) (Mat c.g.) 

(5) 

h = 

flyy ~ lyyG 

V I yyO ) 

) (M at c.g.) 

(6) 

/ _ 

(hz-hzG^ 

) (Mat c.g.) 

(7) 

J 1 ~ 

V hzG J 

h = 

i ^ 

) (M at c.g.) 

(8) 

7 _ 

/u-w 

) (Jyz at C -g-) 

(9) 

J 9 ~ 

V lyzG ) 

7 — 

f hx ~ hxG 

j (U at c.g.) 

(10) 

‘MO — 

V hxG . 


B. Mass Orthogonality Submodule 

This submodule computes each off-diagonal term of the or- 
thogonality matrix using Eq. (11): 

M = Aset Aset (order = nxn) (11) 

In Eq. (11), the order of a reduced-order mass matrix M A , is based on 
the A-set degrees of freedom [ 18 ] and computed using the system 
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equivalent reduction expansion process (SEREP) shown in Eqs. (12) 
and (13) [26]: 


M a = (order = n AD x n AD ) (12) 


= (®Let®Aset) ’^Let ( order = m X n AD ) (13) 


where the matrix 0 Aset (order = n AD X m ) is an eigenmatrix ob- 
tained from Nastran modal analysis, and m and n AD are the number of 
modes to be use in the SEREP and the number of A-set degrees of 
freedom, respectively. The expanded GVT mode shape 0 GAset , is 
also computed using SEREP as follows. 


Reorder the matrix <D Aset , then <I> Aset = 



and the matrix 


0 GAse t (order = n AD X n ) will become Eq. (14): 


^GAset 





(14) 


where matrices (order = n AD X m) and [order = 

(n AD ~ n MD) x m )] are the master and slave eigenmatrices. The 
n MD and n are the measured (or master) degrees of freedom and the 
number of modes for model tuning, respectively; and 0 Gyw (n MD X n) 
is the measured eigenmatrix obtained from GVT. 

The number of performance indices (off-diagonal term square of 
matrix M), 


D. Mode Shape Submodule 

In the previous version of the model-tuning tool, numerical mode 
shapes were improved to minimize mode shape error at sensor 
locations using a single performance index. In the extended version, 
each off-diagonal term of the cross-orthogonality matrix is computed 
using Eq. (20) : 


s = <S> T GAset M A <S> Axt (20) 

The number of performance indices (off-diagonal term square of 
matrix S ) in Eq. (21): 

Jk = sjj i + j i&j = 1,2, ...,n (21) 

computed in this submodule is shown in Eq. (22): 

L s = n 2 — n (22) 

where n is the order of the cross-orthogonality mass matrix S, 
and k = 11 T Lm T tw, . . . , 10 T Lm T m 4- L s • 


E. Modal Assurance Criteria Submodule 

In this submodule, a modal assurance criteria (MAC) matrix 
MA Cij, defined as Eq. (23), 


MACy 


(ELi®L*U 7 


(ELi ®£ g ®hg)(E U kj ) 


(23) 


Jk = M}j i + j i&j = 1,2, (15) 

computed in this submodule is 


is used to compute the corresponding performance indices. The 
diagonal terms of the matrix MACy are selected as m performance 
indices shown in Eq. (24): 


L m — 



(16) 


where k = 11, . . . , 10 + L M . In the previous version of model- 
tuning tool, only one performance index, 


J k = MAC,, (24) 

where i = 1, . . . , m and, k + 11 + L M + m + L s , . . . , 10 + L M + 
2m + L s . 

The optimization problem statement can be written as 


J = E E < 17 > 

i= 1 j=i + 1 

was defined for mass orthogonality. 

The extended version of the model-tuning tool required a 
significant number of performance indices because individual off- 
diagonal terms, as shown in Eqs. (15) and (1A), were performance 
indices instead of using a single norm value of off-diagonal terms as 
in Eq. (17) , which was used in the previous version of the model- 
tuning tool [2,16]. With the extended version, it was much easier to 
improve each off-diagonal term element as opposed to the previous 
version. 

C. Frequency Submodule 

Performance indices associated with the ith frequency error are 
computed from Eq. (18): 

Jk — j- ^ — 11 + L m , . . . , 10 + L m + m (18) 

where/, and f iG are the ith numerical and measured frequencies. The 
m performance indices (number of modes to be matched) are 
computed in this submodule. Similarly, the previous version of the 
model-tuning tool had only one performance index for frequency 
error and was defined as Eq. (19): 



minimize J=^w t Ji (performance index i selected for objective functions) 

such that J k < e k (performance index A; selected for constraint functions) 

(25) 

A constrained optimization problem in Eq. (25) can be solved using 
the O 3 tool. In Eq. (25), e k and w t are small values that can be adjusted 
according to the upper limit values of constraints and the ith 
weighting factor of the performance index / , respectively. 



Fig. 4 X-56A aircraft with soft suspension system. (Photo courtesy 

AFRL/Lockheed Martin). 
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Table 1 Measured and computed flexible modes of the X-56A aircraft 



GVT data 




Nastran results 



Target error, % 

Final design 



Baseline 


Mode number Mode shape 

Frequency 

Frequency 

Error, % 

Mode number 

Frequency 

Error, % 




EFEW configuration 






7 

SW1B 

1.067 

1.035 

-3.0 


7 


1.090 

2.1 

5 

8 

AW1B 

1.543 

1.534 

-0.5 


8 


1.540 

-0.2 

5 

9 

SW1T 

3.223 

2.781 

-13.7 


9 


3.159 

-2.0 

5 

11 

AW1T 

3.839 

3.522 

-8.3 


11 


3.636 

-5.3 

5 

12 

SW2B 

4.440 

4.127 

-7.1 


12 


4.514 

1.7 

5 

13 

AMLGL 

4.466 

4.262 

-4.6 


13 


4.567 

2.3 

5 

14 

SMLGL 

4.666 

4.467 

-4.3 


14 


4.961 

6.3 

5 

15 

BoomH 

5.273 

4.530 

-14.1 


15 


5.223 

-0.9 

10 

18 

AW2B 

6.026 

5.404 

-10.3 


18 


6.061 

0.6 

10 

19 

SWL 

6.264 

5.815 

-7.2 


19 


6.189 

-1.2 

10 

25 

SW3B 

9.346 

9.798 

4.8 


25 


9.416 

0.8 

10 

26 

AW3B 

10.598 

9.889 

-6.7 


27 


11.048 

4.2 

10 

28 

AMLGFA 

11.930 

10.969 

-8.1 


26 


10.035 

-15.9 

10 

30 

AW2T 

12.405 

11.986 

-3.4 


30 


12.811 

3.3 

10 

Total 

weight 

366.7 

— 

— 



366.0 


-0.18 

5 

xc.g. 

location 

165.0 

— 

— 



164.7 


-0.16 

5 

y c.g. 

location 

-0.1 

— 

— 



0.3 


-413 


z C.g. 

location 

N/A 

— 

— 



101.9 


N/A 

— 




FFFW configuration 






7 

SW1B 

1.000 

0.937 

-6.3 


7 


1.001 

0.1 

5 

8 

AW1B 

1.411 

1.392 

-1.3 


8 


1.398 

-0.9 

5 

9 

SW1T 

2.938 

2.608 

-11.2 


9 


2.912 

-0.9 

5 

11 

AW1T 

3.651 

2.932 

-19.7 


11 


3.454 

-5.4 

5 

12 

SW2B 

4.346 

3.898 

-10.3 


12 


4.285 

-1.4 

5 

13 

AMLGL 

4.408 

5.393 

22.4 


13 


4.446 

0.9 

5 

14 

SMLGL 

4.601 

4.159 

-9.6 


14 


4.944 

7.4 

5 

16 

BoomH 

5.276 

4.476 

-15.2 


16 


5.217 

-1.1 

10 

19 

SWL 

6.144 

5.251 

-14.5 


19 


6.018 

-2.0 

10 

24 

SW3B 

8.657 

8.161 

-5.7 


24 


8.673 

0.2 

10 

25 

NLGFA 

9.129 

9.816 

7.5 


25 


9.186 

0.6 

10 

28 

AW2T 

11.540 

10.076 

-12.7 


30 


11.704 

1.4 

10 

Total weight 

488.9 

— 

— 



489.1 


0.04 

5 

v c.g. 

location 

165.2 

— 

— 



165.3 


0.04 

5 

y c.g. 

location 

0.4 

— 

— - — 



0.2 


-41.5 

— 

zc.g. 

location 

N/A 

— 

— 



101.4 


N/A 

— 


III. Ground -Vibration Tests 

The GVT of the X-56A aircraft with a flexible wing configuration 
was completed by ATA Engineering, Inc., under the supervision of 
LMSW, on 17 December 2012. The free-free boundary condition 
with a soft suspension system, as shown in Fig. 4, was used in this 


Table 2 Acronyms for the description of the mode shape 


Acronym 

Mode shape 

AMLGFA 

Antisymmetric main landing gear forward and aft 

AMLGL 

Antisymmetric main landing gear lateral 

AW1B 

Antisymmetric wing first bending 

AW1T 

Antisymmetric wing first torsion 

AW2B 

Antisymmetric wing second bending 

AW2T 

Antisymmetric wing second torsion 

AW3B 

Antisymmetric wing third bending 

AWFA 

Antisymmetric wing forward and aft 

AWL 

Antisymmetric winglet 

BoomH 

Boom horizontal 

BoomV 

Boom vertical 

NLGFA 

Nose landing gear forward and aft 

NLGL 

Nose landing gear lateral 

SMLGFA 

Symmetric main landing gear forward and aft 

SMLGL 

Symmetric main landing gear lateral 

SWL 

Symmetric winglet 

SW1B 

Symmetric wing first bending 

SW1T 

Symmetric wing first torsion 

SW2B 

Symmetric wing second bending 

SW2T 

Symmetric wing second torsion 

SW3B 

Symmetric wing third bending 

SWFA 

Symmetric wing forward and aft 


test. The assembled X-56A aircraft with flexible wings was tested in 
two weight configurations: 1) empty fuel, empty water (EFEW), and 
2) full fuel, full water (FFFW) [1]. It should be noted that, during 
ground vibration and flight testing, water was used to simulate 
fuel weight in the wings. Nondimensional frequency results are 
summarized in Table J_. All the frequencies in this paper were non- 
dimensionalized with respect to the first measured frequency ob- 
tained from GVT with the FFFW configuration. A total of 120 
degrees of freedom have been measured during the GVT. This 
measured mode shape information was expanded to the total of 1477 
A-set degrees of freedom using Eq. (14), and mode shape information 
was computed from Nastran modal analysis of the LMSW’s final 
design model and baseline model with both the EFEW and FFFW 
configurations. A total of the first 30 modes including the six rigid 
body modes were used in this mode shape fitting. In this study, 
the final design model and the baseline model are defined as the 
nonvalidated and the test-validated FE models, respectively. The 
test-validated FE model in Table 1 for the EFEW and the FFFW 
configurations were created by LMSW. Acronyms for the description 
of the mode shape are given in Table 2. The baseline model will be the 
starting configuration of the model-tuning procedure. Figure 5 shows 
how the models are defined in this study. 


IV. Pretuning Analysis 

The GVTs of the X-56A aircraft were performed for the verifi- 
cation of the FE model. The following modal survey requirements are 
observed in [27,28]. 

Military standard : 

1) Analytical model frequencies are to be within 3% of test 
frequencies [27]. 
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$ ■ i 

Fig. 6 Finite-element model of the X-56A aircraft. 



Fig. 5 Different models in this study. NASA standard : 

1) Agreement between test and analysis natural frequencies shall, 
2) Using a cross-orthogonality matrix formed from the analytical as a goal, be within 5% for the significant modes [28]. 
mass matrix and the analytical and test modes, corresponding modes 2) Accurate mass representation of the test article shall be 

are to exhibit at least 95% correlation, and dissimilar modes are to be demonstrated with orthogonality checks using the analytical mass 

orthogonal to within 10% [27]. matrix M and the test mode shapes O g . The orthogonality matrix is 



GVT mode number 

Fig. 7 Comparisons of nondimensional frequencies (EFEW case). 



7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 

GVT mode number 

Fig. 8 Comparisons of nondimensional frequencies (FFFW case). 
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Fig. 9 First five flexible modes with the EFEW configuration. 




QW 1 13 

EFEW cc 

A WT 1 13 

mfiguration 

QW 1 T 

A W 1 T 

OWOT 3 

A A/TT n. T 

ca/tt r:i 

Mode shape 

Mode shape 
Mode number 

oWlt) 

7* 

A W 1 Jj 
8 

oWll 

9 

AW 1 1 
11 

>3 W Zl3 
12 

AlVlLvJL 

13 

ij IVlLvj L 

14 

SW1B 

7 ** 

0.99 

0.00 

0.00 

0.00 

0.08 

0.00 

0.01 

AW1B 

8 

0.00 

1.00 

0.00 

0.00 

0.00 

0.09 

0.00 

SW1T 

9 

0.01 

0.00 

0.91 

0.00 

0.01 

0.00 

0.01 

AW1T 

11 

0.00 

0.01 

0.00 

0.97 

0.00 

0.16 

0.00 

SW2B 

12 

0.09 

0.00 

0.09 

0.01 

0.95 

0.02 

0.27 

AMLGL 

13 

0.00 

0.06 

0.00 

0.43 

0.03 

0.93 

0.02 

SMLGL 

14 

0.05 

0.00 

0.07 

0.00 

0.93 

0.00 

0.51 

FFFW configuration 


Mode shape 

SW1B 

AW1B 

SW1T 

AW IT 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7* 

8 

9 

11 

12 

13 

14 

SW1B 

7 ** 
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0.00 

0.02 

0.00 

0.08 

0.00 

0.01 

AW1B 

8 

0.00 

1.00 

0.00 

0.01 

0.00 

0.07 

0.00 

SW1T 

9 

0.04 

0.00 

0.93 

0.00 

0.01 

0.00 

0.02 

AW1T 

11 

0.00 

0.04 

0.00 

0.90 

0.01 

0.09 

0.01 

SW2B 

12 

0.07 

0.00 

0.06 

0.00 

0.99 

0.00 

0.14 

AMLGL 

13 

0.00 

0.04 

0.00 

0.37 

0.01 

0.96 

0.01 

SMLGL 

14 

0.03 

0.00 

0.04 

0.00 

0.84 

0.02 

0.38 


*: Nastran mode number (Nastran mode extraction order: 1 - 6 = rigid body & higher than 7 = flexible) 
**: GVT mode number (GVT mode extraction order: 1 - 6 = rigid body & higher than 7 = flexible) 


Fig. 10 MAC matrix of the baseline model. 
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1.000 

-0.077 

SMLGL 

14 
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**: GVT mode number 

Fig. 11 Orthogonality matrix of the baseline model. 


computed as O^MO G . As a goal, the off-diagonal terms of the 
orthogonality matrix should be less than 0. 1 for significant modes 
based on the diagonal terms normalized to 1.0 [28]. 

3) Mode shape comparisons shall be required via cross- 
orthogonality checks using the test modes O g , the analytical modes 
O a , and the analytical mass matrix M. The cross-orthogonality 
matrix is computed as O g MO a . As a goal, the absolute value of the 
cross-orthogonality between corresponding test and analytical mode 
shapes should be greater than 0.9, and all other terms of the matrix 
should be less than 0.1 for all significant modes. Additionally, 
qualitative comparisons between test modes and analytical modes 
using mode shape animation and/or deflection plots shall be 
performed [28]. 

Modal analyses of the X-56A aircraft with the EFEW and FFFW 
weight configurations are performed and summarized in the follow- 
ing sections. 

A. Modal Analysis 

For the structural dynamic FE model tuning, the EFEW and FFFW 
weight configurations are selected in this study. A structural dynamic 
FE model of the X-56A aircraft, with a total of 8262 nodes, for 
Nastran modal analysis is shown in Fig. 6. The most important 


frequencies of the X-56A aircraft are summarized in Table 1_ as well 
as Figs. 7, 8. Mode shapes of the first five flexible modes with the 
EFEW configuration are shown in Fig. 9. 

The symmetric and antisymmetric engine lateral modes were not 
captured numerically in this modal analysis with the final design 
model due to the fact that the engines are connected to the center body 
through the use of the rigid bar elements in MSC/Nastran. It should be 
noted that surprisingly large frequency errors are observed with the 
final design model, as shown in Table L Frequency errors for the final 
design model were improved when LMSW performed their model 
tuning. As shown in Table 1, frequency errors of both the EFEW and 
FFFW configurations decreased significantly after LMSW’s model 
tuning, except for a few modes. Frequency errors of mode numbers 
11 and 14 for the EFEW and FFFW configurations and 28 for the 
EFEW configuration are still larger than 5 and 10%, respectively. 

The MAC matrices of the baseline X-56A model with the EFEW 
and FFFW configurations are shown in Fig. 10. In general, diagonal 
terms of the MAC matrices are acceptable because they are above 
90%, with the exception of mode number 14, the symmetric main 
landing gear lateral (SMLGL) mode, for the FFFW configuration. 

The orthogonality and cross-orthogonality matrices of the baseline 
model with the EFEW and FFFW configurations are given in Figs. IT, 
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*: Nastran mode number 
**: GVT mode number 


Fig. 12 Cross-orthogonality matrix of the baseline model (check mode shapes). 
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Fig. 13 Unsteady aerodynamic model of the X-56A aircraft. 


12 . Most of the off-diagonal terms of orthogonality and cross- 
orthogonality matrices were less than 10%; however, a few off-diagonal 
terms were still larger than 10%. From these observations with Table 1 
and Figs. 10-12, it was concluded that LMSW’s baseline model still 
needs to be updated further to have an improvement in accuracy. 

B. Identification of the Primary and Secondary Modes 

One of the primary objectives of the X-56A aircraft is the flight 
demonstration of an active flutter suppression system. In this study, 
the primary and secondary modes for model tuning were identified 
through the use of modal participation factors obtained from a flutter 
analysis of the X-56A aircraft. An unsteady aerodynamic model for 
the flutter analysis using ZAERO code [ 29 ] is shown in Fig. 13. This 
aerodynamic model has 2196 surface elements over the center body, 
the wings, and the winglets. The matched flutter analyses were 
performed at four Mach numbers of 0.130, 0.160, 0.195, and 0.284. 
The speed versus damping V — g and speed versus frequency V - f 
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Fig. 14 V - g and V - f curves of the baseline model with the EFEW at Mach 0.16. 
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curves of the baseline model with the EFEW configuration at Mach 
number of 0.160 from the matched flutter analyses are given in 
Fig. 14. For comparison, the V — g and V - f curves of the final 
design model with the EFEW configuration are also presented in 
Fig. 15. The structural damping used for the flutter speed compu- 
tations was 1%. The flutter speeds and frequencies are summarized in 
Table 3 and Fig. 16. It should be noted in Fig. 16 that the first flutter 
mode is a symmetric body freedom flutter (SBFF) (dashed line), the 
second flutter mode is a symmetric wing bending torsion flutter 
(SWBTF) (dashed single dot line), and the third flutter mode is an 
antisymmetric wing bending torsion flutter (AWBTF) (dashed 
double dot line). 

Flutter mode shapes at Mach numbers of 0.130 and 0.284 are 
shown in Fig. 17. In Fig. 17, center body longitudinal motions 
together with the outboard wing bending motion can be observed in 
the first SBFF motion. These outboard wing bending motions are 


more of a divergence type of motion or wash in motion. On the other 
hand, small center body pitch and wing bending in addition to torsion 
motions are noticed in the first symmetric flutter mode. Small center 
body roll and outboard wing wash out bending and torsion motion 
can be seen in the case of the first antisymmetric flutter mode. An 
antisymmetric body freedom flutter mode is captured with a Mach 
number of 0.284 (open diamond marker in Fig. 16), and the 
corresponding mode shape is also given in Fig. 17 . 

Modal participation factors are summarized in Table 4. In the case 
of the final design model, modal participation of mode numbers 7,8, 
9, and 1 1 together with rigid body modes are 89.0-96.5% for all three 
Mach number cases as shown in Table 4. Therefore, these modes 
were the primary modes for the final design model. However, in the 
case of the baseline model, modal participation of the same modes 
was reduced to 84.7-93.0% for the SBFF mode, 76.3-89.5% for the 
symmetric flutter mode, and 33.6-68.3% for the antisymmetric 




Fig. 15 V — g and V — f curves of the final design model with the EFEW at Mach 0.16. 
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Final design (non-validated) 

• Symmetric body freedom flutter 

■ Symmetric wing bending torsion flutter 

♦ Anti-symmetric wing bending torsion flutter 


Baseline (validated by LMSW) 

O Symmetric body freedom flutter= 

□ Symmetric wing bending torsion flutter 
O Anti-symmetric wing bending torsion flutter 


After model tuning (validated by NASA) 

© Symmetric body freedom flutter — > Model tuning effect (by LMSW) 

0 Symmetric wing bending torsion flutter — > Model tuning effect (by NASA) 



Fig. 16 Flutter boundaries and flight envelope under the EFEW configuration. 


flutter mode. Therefore, mode numbers 12, 13, and 14 for the EFEW 
configuration and the same mode numbers for the FFFW configura- 
tion were added to the previous primary mode to have modal 
participation of 93 .3-98.9%. It should be noted that mode number 14, 
the symmetric main landing gear lateral mode, involves an average of 
7.4 and 8.4% for the symmetric wing bending torsion flutter in the 
case of the EFEW, and for the antisymmetric wing bending torsion 
flutter in the case of the FFFW, respectively. Other than these two 
cases, this mode contributes less than 2.6% of the symmetric body 
freedom flutter mode. 


V. Structural Dynamic Model Tuning 

In general, the major objective of a GVT is the modal validation of 
the structural dynamic FE model. Based on the results of the test and 
analysis correlation, if these correlation results violate the military 
and NASA standards [27,28], then the FE model needs to be adjusted 
to match the test data. 

In this study, based on the data shown in Table _l and Figs. 10-12, 
the baseline model developed by LMSW is further updated to 
have improved correlation with test data. The four model-tuning 
procedures performed in this study used the DOT algorithm to save 


Antisymmetric flutter at Mach 0.130 
• Small center body roll+ wing bending and torsion 
(wash out motion) 



Fig. 17 First three flutter mode shapes of the X-56A baseline model. 
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T able 3 Nondimensional flutter speeds and frequencies with EFEW 


at Mach 0.16 

Baseline 

After fourth tuning 

Flutter mode Final design Speed Difference, 2 % 

Speed Difference^ % 


Nondimensional flutter speed 


First 

0.981 

1.135 

15.7 

1.190 

21.3 

Second 

1.090 

1.462 

34.1 

1.580 

44.9 

Third 

1.161 

1.809 

55.8 

1.883 

62.2 


Nondimensional flutter frequency 


First 

0.653 

0.663 

1.6 

0.664 

1.7 

Second 

2.014 

2.415 

19.9 

2.527 

25.5 

Third 

1.390 

2.509 

80.5 

2.755 

98.2 


difference = (baseline-final design) /final design. 
b Difference = (after tuning-final design) /final design. 

computation time. A summary of target objective functions to be 
improved for each model-tuning procedure are shown in Fig. 1_8. 

A. First Model-Tuning Procedure 

Frequency errors of the primary mode numbers 1 1 and 14 for the 
EFEW and FFFW configurations, with frequency errors of —5.3, 6.3, 
-5.4, and 7.4% in Table 1, are selected as an objective function. The 
two modes, mode number 14 for the EFEW and FFFW configura- 
tions, are mainly related to the main landing gear. Therefore, one 
lumped mass value and nine sectional properties of the main landing 
gear, together with the corresponding Young’s modulus E and shear 
modulus G were selected as design variables. It was assumed to have 
symmetric structural properties, and therefore design variable linking 
was used for these design variables. In this study, lumped mass values 
of the GVT sensor cables, 8 lb in total, are also selected as design 
variables. 

Differences in the total weights (1 + 1), x c.g. locations (1 + 1), all 
other frequencies (13 + 10), the off-diagonal terms of symmetric 
orthogonality matrices (21 +21) in Eq. (13), and the off-diagonal 
terms of cross-orthogonality matrices (42 + 42) in Eq. (14) were 
used as inequality constraints. Numbers in parenthesis designate that 
the first and second numbers are the number of inequality constraints 
obtained from the EFEW and FFFW configurations, respectively. 
Therefore, the total number of 12 design variables as well as 153 
inequality constraints (2 + 2 + 23 + 42 + 84) were used in this first 
optimization procedure. In this study, MAC values were not used for 
the constraint functions. Plus and minus 20% of the starting design 
variable values are selected as the upper and lower bounds for design 
variables (i.e., side constraints), respectively. 

Continuous design variables were assumed, and model tuning was 
based on the DOT algorithm in the O 3 tool. Updated frequency results 
and the MAC matrix are shown in Table 5 and Fig. 19, respectively. 
Orthogonality and cross-orthogonality matrices for the EFEW and 
FFFW configurations after the first model-tuning procedure are given 
in Figs. 20, 21 . After the first model tuning, four frequencies related 
to the objective function improved from -5.3, 6.3, -5.4, and 7.4% to 
-3.5, 5.6, -3.9, and 6.7%, respectively. Mode numbers 1 1 and 14 for 
the EFEW and FFFW configurations after the first model tuning still 
violate the military standard for primary modes (3%). However, 
mode number 1 1 for these two weight configurations satisfied the 
NASA standard for primary modes (5%). 

In Fig. 20, small improvements for the off-diagonal terms of the 
orthogonality matrix were observed; 0.143 and 0.150 became 0.140 
and 0.146, respectively. On the other hand, -0.143 for the EFEW 
case became -0.144, which was mainly because, when upper limit 
values of the constraint were prepared, these values were rounded up 
-0.143 to -0.144 to give some buffer value. This off-diagonal term 
was an active constraint, and therefore the optimizer tried to increase 
this constraint function during the first optimization procedure. 

Big improvements in off-diagonal terms of the cross-orthogonality 
matrix were observed as shown in Fig. 21. In addition, a few off- 
diagonal terms were approaching the 10% target value. It should be 
noted that four off-diagonal terms violated military and NASA 
standards; -24.4, 13.3,-13.3, and -10.9%, in Fig. 12, were less than 


10% target value after the first model tuning, and 9.2, 2.6, -11.0, and 
-7.6%, in Fig. 21. 

It should be emphasized that the frequencies of the mode numbers 
12 and 13 for the EFEW weight configuration were too close, as 
shown in Table 5, and the frequency difference was less than 0.001 
(nondimensional frequency). As shown in Fig. 19, these two modes 
were not interchanged yet. Once mode interchange happens during 
a model-tuning procedure, then all the off-diagonal terms of the 
orthogonality and cross-orthogonality matrices will be discontinu- 
ous, and therefore the sensitivity matrices computed for finding the 
search direction and surface slope changes will be corrupted. To 
overcome this mode interchange problem, a special mode-tracking 
program based on diagonal and off-diagonal terms of the MAC 
matrix was developed and implemented into the current extended 
model-tuning tool. This program helped with the continuity of the 
performance indices. 

B. Second Model-Tuning Procedure 

For the second tuning procedure, mode numbers 12 and 14 of the 
cross-orthogonality matrix for the EFEW weight configuration 
were selected as the objective functions; specifically, the objective 
functions selected were off-diagonal terms 12-14 and 14-12. From 
Fig. 21, after the first model-tuning procedure, these off-diagonal 
terms resulted in an error of nearly 36 and 50%, respectively, 
and needed to be improved further. It should be noted that four 
improved frequencies, mode numbers 1 1 and 14 for the two weight 
configurations, were switched to the constraint functions for this run. 
In this procedure, 20 design variables (eight from the last procedure, 
four from GYT sensor cable masses, and eight additional from main 
landing gear properties) were selected from a sensitivity analysis. 
This decision making is due to the fact that mode number 14 for both 
weight configurations is related directly to the main landing gear. In 
addition, upper and lower bounds for these design variables were kept 
at ±20%. All other cross-orthogonality matrix terms (40 + 42) were 
used as constraint functions along with the remainder of the 
orthogonality matrix terms (21 +21), frequency terms (15 + 12), 
total weight (1 + 1), x c.g. locations (1 + 1), and GVT sensor cable 
weight (1) for a total of 156 constraint functions. Like the previous 
run, no MAC constraints were used. 

It should be noted here that mode number 12 of the EFEW weight 
configuration can be interchanged with mode number 13 during the 
second model-tuning procedure. Actually, the main reason why a 
special mode-tracking routine was developed in this study was due to 



Fig. 18 Summary of objective functions. 
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Table 4 Modal participation factors (percent) (first, second, and third flutter modes) 


Final design 


Baseline model 


GVT mode number 

Mode shape 


M = 0.130 


M = 0.195 

M = 0.284 


M = 0.130 


M = 0.195 

M = 0.284 

First 

Second 

Third 

First 

Second 

Third 

First 

Second 

Third 

First 

Second 

Third 

First 

Second 

Third 

First 

Second 

Third 









EFEW configuration 











Primary modes 

1-6 

Rigid 

31.6 

30.7 

40.3 

27.5 

32.6 

33.0 

25.0 

34.3 

25.2 

33.7 

33.3 

42.8 

27.9 

35.7 

40.5 

24.6 

39.0 

40.2 


7 

SW1B 

15.0 

9.5 

0.0 

12.1 

8.8 

0.0 

9.7 

8.1 

0.0 

17.0 

10.0 

0.0 

14.9 

9.2 

0.0 

13.0 

8.6 

0.0 


8 

AW1B 

0.0 

0.0 

27.3 

0.0 

0.0 

31.1 

0.0 

0.0 

35.1 

0.0 

0.0 

8.3 

0.0 

0.0 

12.5 

0.0 

0.0 

28.1 


9 

SW1T 

44.3 

54.6 

0.0 

51.1 

54.4 

0.0 

56.1 

54.1 

0.0 

38.6 

43.0 

0.0 

47.8 

41.5 

0.0 

53.7 

39.7 

0.0 


11 

AW1T 

0.0 

0.0 

27.3 

0.0 

0.0 

31.1 

0.0 

0.0 

35.1 

1.9 

2.8 

0.0 

1.8 

2.5 

0.0 

1.7 

2.2 

0.0 


12 

SW2B 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

42.1 

0.0 

0.0 

40.9 

0.0 

0.0 

28.1 


13 

AMLGL 

1.5 

1.3 

0.0 

1.3 

0.9 

0.0 

1.2 

0.7 

0.0 

0.0 

0.0 

4.5 

0.0 

0.0 

4.2 

0.0 

0.0 

2.5 


14 

SMLGL 

1.3 

0.7 

0.0 

1.2 

0.6 

0.0 

1.2 

0.6 

0.0 

2.6 

7.1 

0.0 

1.7 

7.6 

0.0 

1.2 

7.5 

0.0 



Total 

93.7 

96.8 

95.0 

93.2 

97.3 

95.2 

93.2 

97.8 

95.4 

93.8 

96.2 

97.7 

94.1 

96.5 

98.1 

94.2 

97.0 

98.9 

Secondary modes 

15 

BoomH 

0.0 

0.0 

0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.7 

1.6 

0.0 

0.7 

1.7 

0.0 

0.6 

1.7 

0.0 


26 

AMLGFA 

0.0 

0.0 

1.6 

0.0 

0.0 

1.6 

0.0 

0.0 

1.7 

0.0 

0.0 

0.4 

0.0 

0.0 

0.4 

0.0 

0.0 

0.4 


28 

SW2T 

2.9 

1.1 

0.0 

3.1 

0.9 

0.0 

3.2 

0.8 

0.0 

2.6 

1.0 

0.0 

2.5 

0.8 

0.0 

2.5 

0.6 

0.0 


30 

AW2T 

0.0 

0.0 

1.6 

0.0 

0.0 

1.9 

0.0 

0.0 

1.9 

1.8 

0.5 

0.0 

1.7 

0.4 

0.0 

1.7 

0.3 

0.0 



Total 

2.9 

1.1 

3.3 

3.1 

0.9 

3.5 

3.2 

0.8 

3.6 

4.1 

3.1 

1.4 

4.9 

2.9 

0.4 

4.8 

2.6 

0.4 









FFFW configuration 











Primary modes 

1-6 

Rigid 

42.4 

32.3 

44.7 

38.5 

36.2 

38.4 

34.8 

40.5 

39.0 

42.1 

29.2 

35.3 

36.8 

32.0 

34.4 

32.6 

36.0 

32.4 


7 

SW1B 

12.9 

11.5 

0.0 

11.8 

10.9 

0.0 

11.1 

10.3 

0.0 

14.9 

10.4 

0.0 

12.8 

9.5 

0.0 

10.7 

8.9 

0.0 


8 

AW1B 

0.0 

0.0 

5.2 

0.0 

0.0 

27.2 

0.0 

0.0 

25.0 

0.0 

0.0 

1.6 

0.0 

0.0 

1.3 

0.0 

0.0 

1.2 


9 

SW1T 

38.0 

46.3 

0.0 

42.0 

42.9 

0.0 

45.9 

39.5 

0.0 

29.7 

37.1 

0.0 

35.4 

34.1 

0.0 

40.6 

30.8 

0.0 


11 

AW1T 

0.0 

0.0 

44.0 

0.0 

0.0 

27.2 

0.0 

0.0 

25.0 

0.7 

0.8 

0.0 

0.7 

0.7 

0.0 

0.8 

0.7 

0.0 


12 

SW2B 

0.1 

0.1 

0.0 

0.1 

0.1 

0.0 

0.1 

0.1 

0.0 

0.0 

0.0 

50.2 

0.0 

0.0 

50.1 

0.0 

0.0 

50.7 


13 

AMLGL 

0.0 

0.0 

2.5 

0.0 

0.0 

1.3 

0.0 

0.0 

1.4 

7.6 

20.3 

0.0 

8.2 

21.9 

0.0 

8.6 

22.2 

0.0 


14 

SMLGL 

1.9 

6.8 

0.0 

1.6 

7.4 

0.0 

1.3 

7.5 

0.0 

0.0 

0.0 

7.3 

0.0 

0.0 

8.4 

0.0 

0.0 

9.5 



Total 

95.3 

97.0 

96.4 

94.0 

97.5 

94.1 

93.2 

97.9 

90.4 

95.0 

97.8 

94.4 

93.9 

98.2 

94.2 

93.3 

98.6 

93.8 

Secondary modes 

16 

BoomH 

0.0 

0.0 

1.1 

0.0 

0.0 

1.2 

0.0 

0.0 

1.4 

0.0 

0.0 

1.8 

0.0 

0.0 

1.9 

0.0 

0.0 

2.1 


19 

SWL 

0.5 

0.7 

0.0 

0.7 

0.6 

0.0 

0.8 

0.5 

0.0 

0.0 

0.0 

2.3 

0.0 

0.0 

2.3 

0.0 

0.0 

2.4 


24 

SW3B 

1.3 

0.5 

0.0 

1.6 

0.5 

0.0 

1.9 

0.5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


25 

NLGFA 

2.0 

0.9 

0.0 

2.5 

0.7 

0.0 

2.9 

0.6 

0.0 

3.1 

1.2 

0.0 

3.8 

0.9 

0.0 

4.3 

0.8 

0.0 


30 

AW2T 

0.0 

0.0 

1.1 

0.0 

0.0 

3.0 

0.0 

0.0 

5.6 

0.9 

0.2 

0.0 

1.1 

0.2 

0.0 

1.2 

0.2 

0.0 



Total 

3.8 

2.1 

2.2 

4.8 

1.8 

4.2 

5.6 

1.6 

7.0 

4.0 

1.4 

4.1 

4.9 

1.1 

4.2 

5.5 

1.0 

4.5 
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EFEW configuration 


Mode shape 

SW1B 

AW1B 

SW1T 

AW1T 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7* 

8 

9 

11 

12 

13 

14 

SW1B 

7** 

0.99 

0.00 

0.00 

0.00 

0.08 

0.00 

0.01 

AW1B 

8 

0.00 

1.00 

0.00 

0.00 

0.00 

0.09 

0.00 

SW1T 

9 

0.01 

0.00 

0.92 

0.00 

0.00 

0.00 

0.01 

AW1T 

11 

0.00 

0.01 

0.00 

0.98 

0.01 

0.17 

0.00 

SW2B 

12 

0.09 

0.00 

0.07 

0.01 

0.98 

0.00 

0.30 

AMLGL 

13 

0.00 

0.06 

0.00 

0.41 

0.00 

0.97 

0.02 

SMLGL 

14 

0.05 

0.00 

0.06 

0.00 

0.92 

0.03 

0.55 



cvxr 1 o 

FFFW cc 

A W 7 1 R 

infiguration 
cvxn "T 

A W 1 T 

QYX/OR 

A A/TT m 

c A /TT rZT 

Mode shape 

Mode shape 
Mode number 

o W lr> 

7* 

AW IrS 
8 

5 W 1 1 

9 

AW I 1 
11 

OWZD 

12 

AlVILvj L 

13 

oiVILLiL 

14 

SW1B 

7** 

0.99 

0.00 

0.02 

0.00 

0.08 

0.00 

0.01 

AW1B 

8 

0.00 

1.00 

0.00 

0.02 

0.00 

0.07 

0.00 

SW1T 

9 

0.04 

0.00 

0.95 

0.00 

0.00 

0.00 

0.02 

AW1T 

11 

0.00 

0.04 

0.00 

0.93 

0.01 

0.10 

0.01 

SW2B 

12 

0.07 

0.00 

0.04 

0.00 

0.99 

0.00 

0.15 

AMLGL 

13 

0.00 

0.04 

0.00 

0.33 

0.01 

0.97 

0.01 

SMLGL 

14 

0.03 

0.00 

0.03 

0.00 

0.85 

0.02 

0.38 


*: Nastran mode number 
**: GVT mode number 


Fig. 19 MAC matrix after the first model-tuning procedure (DOT-01). 


this mode interchange. The MAC matrix after the second model- 
tuning procedure is given in Fig. 22. 

After several optimization iterations, it turned out that the fre- 
quency errors of mode number 14 (SMLGL) and mode number 28 
(AMLGFA) for the EFEW configuration, off-diagonal term 12-13 
(AMLGL SW2B) of the orthogonality matrix for the EFEW con- 
figuration, and finally off-diagonal term 14-13 (AMLGL-SMLGL) 
of the cross-orthogonality matrix became active constraints, and 
no improvements were observed. Therefore, the constraint limit 
values were relaxed to give more allowance for improvement, 6 
and 20% upper limit for frequency errors, 15% for off-diagonal terms 
12-13 and 14-13 of the orthogonality and cross-orthogonality 
matrices. 

Model tuning was performed using the DOT algorithm. Frequency 
results for this second tuning procedure are also shown in Table 5. 
Orthogonality and cross-orthogonality matrices for the EFEW 
and FFFW configurations are presented in Figs. 23, 24. For both 
configurations, frequency errors did not change significantly. Even 
with relaxed upper limit values, frequency errors for modes 14 and 28 
did not change at all. 


Results for the orthogonality matrices for both weight configu- 
rations after the second procedure are shown in Fig. 23 . Note that, for 
this procedure, all of the off-diagonal terms for the orthogonality 
matrices were selected as constraints, not objective functions. 
Similarly, orthogonality matrices for both weight configurations 
were unchanged during the second optimization run. Off-diagonal 
term 12-13 was an active constraint, and this term did not become 
better or worse even with relaxed upper limit value. 

Results for the off-diagonal terms of the cross-orthogonality 
matrix after the second tuning procedure are shown in Fig. 24. In 
comparison to the first tuning procedure results in Fig. 21, some 
off-diagonal terms, mostly constraints, slightly improved for both 
configurations. However, the objective functions for this procedure 
were off-diagonal terms 12-14 and 14-12 for the EFEW configu- 
ration and have improved negligibly from -0.497 to -0.495 and 
from -0.358 to -0.357. 

C. Third Model-Tuning Procedure 

Based on the results from the second model-tuning procedure, 
there was very minimal improvement in the cross-orthogonality 




CW1 R 

EFEW cc 

A W 1 R 

mfiguration 

Q W 1 T 

A W 1 T 

CW 9 R 

aa/TT m 

ca/ft nr 

Mode shape 

Mode shape 
Mode number 

jW Id 
7 ** 

AW 1 D 
8 

jWI I 

9 

AW 1 1 
11 

o W Z o 
12 

AlVlLuL 

13 

1 V 1 J — < J—/ 

14 

SW1B 

7 ** 

1.000 

-0.021 

-0.053 

-0.011 

0.026 

0.024 

-0.033 

AW1B 

8 

-0.021 

1.000 

0.012 

0.002 

0.004 

0.039 

0.002 

SW1T 

9 

-0.053 

0.012 

1.000 

0.004 

0.035 

-0.007 

-0.025 

AW1T 

11 

-0.011 

0.002 

0.004 

1.000 

0.021 

-0.091 

0.003 

SW2B 

12 

0.026 

0.004 

0.035 

0.021 

1.000 

-0.144 

0.140 

AMLGL 

13 

0.024 

0.039 

-0.007 

-0.091 

-0.144 

1.000 

0.005 

SMLGL 

14 

-0.033 

0.002 

-0.025 

0.003 

0.140 

0.005 

1.000 



cwi R 

FFFW co 

A W 1 R 

mfiguration 

Q W 1 T 

A W 1 T 

QW 9 R 

AA/TT fTT 

CA/fT rij 

Mode shape 

Mode shape 
Mode number 

oWlD 

7 ** 

AW 1 D 
8 

oWll 

9 

AW 1 1 
11 

o W ZrS 
12 

AIVILCjL 

13 

uIVILvjL 

14 

SW1B 

7** 

1.000 

0.013 

-0.047 

0.014 

0.019 

0.008 

-0.065 

AW1B 

8 

0.013 

1.000 

-0.010 

0.013 

-0.005 

-0.062 

-0.011 

SW1T 

9 

-0.047 

-0.010 

1.000 

-0.019 

0.007 

-0.013 

-0.026 

AW1T 

11 

0.014 

0.013 

-0.019 

1.000 

-0.026 

0.092 

0.017 

SW2B 

12 

0.019 

-0.005 

0.007 

-0.026 

1.000 

0.003 

0.146 

AMLGL 

13 

0.008 

-0.062 

-0.013 

0.092 

0.003 

1.000 

-0.076 

SMLGL 

14 

-0.065 

-0.011 

-0.026 

0.017 

0.146 

-0.076 

1.000 


**: GVT mode number 

Fig. 20 Orthogonality matrix after the first model-tuning procedure (DOT-01). 
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EFEW configuration 


Mode shape 

SW1B 

AW1B 

SW1T 

AW1T 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7* 

8 

9 

11 

12 

13 

14 

SW1B 

7** 

1.006 

0.015 

0.034 

0.007 

-0.022 

-0.015 

-0.064 

AW1B 

8 

-0.006 

-1.018 

-0.004 

0.008 

0.000 

-0.030 

0.003 

SW1T 

9 

-0.028 

-0.007 

-1.038 

-0.020 

0.106 

0.021 

-0.021 

AW1T 

11 

0.000 

-0.011 

-0.003 

-1.032 

0.001 

-0.064 

-0.002 

SW2B 

12 

0.000 

-0.004 

-0.126 

-0.019 

-0.894 

0.092 

-0.358 

AMLGL 

13 

0.009 

-0.015 

-0.016 

0.180 

0.026 

-0.983 

-0.002 

SMLGL 

14 

0.010 

-0.003 

-0.071 

0.006 

-0.497 

-0.009 

0.759 



O w 1 R 

FFFW cc 

A W 1 R 

infiguration 

cw 1 T 

A W 1 T 

Q W9R 

A A /TT riT 

ca/tt m 

Mode shape 

Mode shape 
Mode number 

o W ir> 

7* 

i\ W IJj 
8 

jW 1 I 

9 

/\W I 1 
11 

o W ZJj 
12 

AiVILvJ L 

13 

ulVIL/vJE/ 

14 

SW1B 

7** 

-1.086 

-0.008 

-0.037 

0.012 

0.033 

0.012 

0.061 

AW1B 

8 

-0.003 

-1.109 

-0.008 

0.011 

-0.004 

-0.059 

0.011 

SW1T 

9 

0.026 

0.004 

1.146 

-0.016 

-0.110 

-0.007 

0.016 

AW1T 

11 

0.019 

0.001 

-0.032 

1.082 

0.014 

-0.076 

-0.009 

SW2B 

12 

0.009 

0.003 

0.112 

-0.005 

1.046 

-0.053 

0.205 

AMLGL 

13 

0.003 

0.010 

0.002 

0.180 

0.072 

1.036 

0.002 

SMLGL 

14 

0.011 

0.011 

0.038 

0.014 

0.340 

-0.113 

-0.843 


*: Nastran mode number 
**: GVT mode number 


Fig. 21 Cross-orthogonality matrix after the first model-tuning procedure (DOT-01). 


matrix as shown in Fig. 24, when compared to the first model-tuning 
run in Fig. 21. The objective functions and corresponding design 
variables during the previous model-tuning procedure were reviewed 
before conducting this third model-tuning procedure. By looking at 
some of the previous designs, there was a possibility that a better 
design configuration existed. Just because the output of the previous 
procedure resulted in miniscule improvement does not imply that the 
O 3 tool program did not find a better optimization configuration in 
one of its iterations. From the previous optimization history (POH), a 
better starting configuration was selected. This design configuration 
was not selected by the DOT because some of the performance 
indices had violated constraint functions, which were actually 
acceptable violations. The results from this selected configuration are 
displayed in Table 6 and Figs. 25-27 for frequencies, the MAC 
matrix, the mass orthogonality matrix, and the cross-orthogonality 
matrix, respectively. 


When comparing frequency results in Tables 5 and 6, it can be seen 
that there was a mode interchange between modes 12 and 13 for both 
weight configurations as the frequencies were switched. As shown in 
Tables 5 and 6, frequencies of mode number 14 for the EFEW and 
FFFW configurations after DOT-02 in Table 5,5.6 and 6.7%, were all 
less than 5% (even smaller than 3%) in Table 6 (under POH). Also, 
notice that the frequency for mode 28 of the EFEW weight 
configuration reached the constraint value limit of 20%. In Fig. 25, 
the MAC values for mode 14 of the EFEW and FFFW configurations 
based on the selected POH was 0.97 and 0.61, respectively; these 
MAC values are far better than the 0.55 and 0.39 from DOT-02 as 
shown in Fig. 22. In addition, the mass orthogonality matrices for 
both weight configurations improved as well, as shown in Fig. 26 ; 
however, for the EFEW, off-diagonal term 12-13 increased when 
compared to the second model-tuning results, whereas off-diagonal 
term 12-14 improved to 0. 107 when compared to 0. 140 in Fig. 23. For 


EFEW configuration 


Mode shape 

SW1B 

AW1B 

SW1T 

AW1T 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7* 

8 

9 

11 

12 

13 

14 

SW1B 

7 ** 

0.99 

0.00 

0.00 

0.00 

0.08 

0.01 

0.01 

AW1B 

8 

0.00 

1.00 

0.00 

0.00 

0.00 

0.08 

0.00 

SW1T 

9 

0.01 

0.00 

0.92 

0.00 

0.00 

0.00 

0.01 

AW1T 

11 

0.00 

0.01 

0.00 

0.98 

0.02 

0.16 

0.00 

SW2B 

12 

0.09 

0.00 

0.07 

0.01 

0.98 

0.01 

0.30 

AMLGL 

13 

0.00 

0.06 

0.00 

0.41 

0.00 

0.95 

0.02 

SMLGL 

14 

0.05 

0.00 

0.06 

0.00 

0.91 

0.07 

0.55 

FFFW configuration 

— 

Mode shape 

SW1B 

AW1B 

SW1T 

AW IT 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7* 

8 

9 

11 

12 

13 

14 

SW1B 

7 ** 

0.99 

0.00 

0.02 

0.00 

0.08 

0.00 

0.01 

AW1B 

8 

0.00 

1.00 

0.00 

0.02 

0.00 

0.07 

0.00 

SW1T 

9 

0.04 

0.00 

0.95 

0.00 

0.00 

0.00 

0.02 

AW1T 

11 

0.00 

0.04 

0.00 

0.93 

0.01 

0.10 

0.01 

SW2B 

12 

0.07 

0.00 

0.04 

0.00 

0.99 

0.00 

0.15 

AMLGL 

13 

0.00 

0.04 

0.00 

0.33 

0.01 

0.97 

0.01 

SMLGL 

14 

0.03 

0.00 

0.03 

0.00 

0.85 

0.02 

0.39 


*: Nastran mode number 
**: GVT mode number 


Fig. 22 MAC matrix after the second model-tuning procedure (DOT-02). 
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EFEW configuration 


Mode shape 

SW1B 

AW1B 

SW1T 

AW1T 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7 ** 

8 

9 

11 

12 

13 

14 

SW1B 

7 ** 

1.000 

-0.021 

-0.053 

-0.011 

0.026 

0.024 

-0.033 

AW1B 

8 

-0.021 

1.000 

0.012 

0.002 

0.004 

0.039 

0.002 

SW1T 

9 

-0.053 

0.012 

1.000 

0.004 

0.035 

-0.007 

-0.025 

AW1T 

11 

-0.011 

0.002 

0.004 

1.000 

0.021 

-0.091 

0.003 

SW2B 

12 

0.026 

0.004 

0.035 

0.021 

1.000 

-0.144 

0.140 

AMLGL 

13 

0.024 

0.039 

-0.007 

-0.091 

-0.144 

1.000 

0.005 

SMLGL 

14 

-0.033 

0.002 

-0.025 

0.003 

0.140 

0.005 

1.000 

FFFW configuration 

— — 

Mode shape 

SW1B 

AW1B 

SW1T 

AW1T 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7 ** 

8 

9 

11 

12 

13 

14 

SW1B 

7 ** 

1.000 

0.013 

-0.047 

0.014 

0.019 

0.008 

-0.065 

AW1B 

8 

0.013 

1.000 

-0.010 

0.013 

-0.005 

-0.062 

-0.011 

SW1T 

9 

-0.047 

-0.010 

1.000 

-0.019 

0.007 

-0.013 

-0.026 

AW1T 

11 

0.014 

0.013 

-0.019 

1.000 

-0.026 

0.092 

0.017 

SW2B 

12 

0.019 

-0.005 

0.007 

-0.026 

1.000 

0.003 

0.146 

AMLGL 

13 

0.008 

-0.062 

-0.013 

0.092 

0.003 

1.000 

-0.076 

SMLGL 

14 

-0.065 

-0.011 

-0.026 

0.017 

0.146 

-0.076 

1.000 


**: GVT mode number 

Fig. 23 Orthogonality matrix after the second model-tuning procedure (DOT-02). 


the FFFW, off-diagonal term 12-14 also significantly improved from 
0.146 to 0.118. 

Furthermore, Fig. 27 displays a strong improvement in the cross- 
orthogonality matrices when compared to Fig. 24 after the second 
model-tuning procedure. After the first model-tuning procedure in 
Fig. 21, off-diagonal terms 12-13 and 13-12 for the EFEW 
configuration were constrained to be a maximum value of 10%. After 
the second model-tuning procedure, in Fig. 24, these two off- 
diagonal terms satisfied the 10% requirement. Although these two 
off-diagonal terms in Fig. 27 violated the 10% requirement; however, 
compared to Fig. 24 results, the cross-orthogonality matrix in Fig. 27 
had a more improved configuration when compared to the results in 
Fig. 24. Significant improvements were observed for the EFEW off- 
diagonal terms 12-14 and 14-12, which went from -0.357 to -0.139 
and -0.495 to 0.000, respectively. Likewise, the FFFW off-diagonal 
terms 12-14 and 14-12 (14-13 due to mode interchange) improved 
magnitudewise from 0.205 to 0. 162 and from 0.340 to 0.286. On the 
other hand, some mode shape off-diagonal terms did not improve; 


most significant was the EFEW off-diagonal term 12-13 (12-12 due 
to mode interchange), which was 0.030 after the second model- 
tuning procedure but was 0.238 based on POH. Overall, this POH 
configuration is much better than the results after the second model- 
tuning procedure. By selecting an improved optimization configu- 
ration based on the POH as the new starting point for this run, the 
optimization results for this third model-tuning procedure can be 
improved significantly. 

Based on the cross-orthogonality matrix from the POH in 
Fig. 27, off-diagonal terms 12-12 (SW2B-AMLGL) of the EFEW 
configuration and 14-13 (SMLGL-SW2B) of the FFFW configu- 
ration were selected as the objective functions for this third model- 
tuning procedure. The same design variables were used as before, 
but with the addition of four main landing gear design variables, 
for a total of 24 design variables. Constraint functions consisted 
of all other cross-orthogonality matrix terms (41 +41), mass 
orthogonality matrix terms (21 +21), frequency terms (15 + 12), 
total weight (1 + 1), v c.g. locations (1 + 1), and GVT sensor cable 


EFEW configuration 

— — 

Mode shape 

SW1B 

AW1B 

SW1T 

AW1T 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7* 

8 

9 

11 

12 

13 

14 

SW1B 

7** 

-1.006 

0.015 

0.034 

0.007 

- 0.022 

-0.017 

-0.064 

AW1B 

8 

0.006 

-1.018 

-0.004 

0.008 

0.002 

-0.030 

0.003 

SW1T 

9 

0.028 

-0.007 

-1.038 

- 0.020 

0.104 

0.028 

- 0.021 

AW1T 

11 

0.000 

- 0.011 

-0.003 

-1.032 

0.005 

-0.064 

- 0.002 
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Fig. 24 Cross-orthogonality matrix after the second model-tuning procedure (DOT-02). 
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Table 5 Flexible modes after the first and second model-tuning procedures 
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weight (_L) for a total of 156 constraint functions. Again, the DOT 
algorithm was used for this optimization procedure. For both weight 
configurations, constraint value limits of 15% were used for 
off-diagonal terms 12-13, 12-14, 13-12, and 14-12 of the 


orthogonality matrices and off-diagonal terms 12-12 (mode 
interchange), 12-14, 13-13 (mode interchange), and 14-13 of the 
cross-orthogonality matrices. For all other off-diagonal terms, 10% 
limit values were used. 




QYX/1 T3 

EFEW cc 
a vxn r> 

mfiguration 

QYX7U T 

A \\T 1 T 


a A/TT nr 

qa/tt n r 

Mode shape 

Mode shape 
Mode number 

oW IrS 

7* 

AW Its 
8 

o W 1 1 

9 

AW 1 I 
11 

15 wzts 

13 

AlVlLCrL 

12 

uIVILvjL 

14 

SW1B 

7** 

0.99 

0.00 

0.00 

0.00 

0.08 

0.00 

0.07 

AW1B 

8 

0.00 

1.00 

0.00 

0.00 

0.00 

0.10 

0.00 

SW1T 

9 

0.01 

0.00 

0.92 

0.00 

0.00 

0.00 

0.01 

AW1T 

11 

0.00 

0.01 

0.00 

0.98 

0.00 

0.29 

0.00 

SW2B 

12 

0.09 

0.00 

0.07 

0.01 

0.95 

0.03 

0.89 

AMLGL 

13 

0.00 

0.06 

0.00 

0.38 

0.03 

0.96 

0.03 

SMLGL 

14 

0.05 

0.00 

0.06 

0.00 

0.88 

0.00 

0.97 

FFFW configuration 

— — 

Mode shape 

SW1B 

AW1B 

SW1T 

AW IT 

SW2B 

AMLGL 

SMLGL 

Mode shape 

Mode number 

7* 

8 

9 

11 

13 

12 

14 

SW1B 

7** 

0.99 

0.00 

0.02 

0.00 

0.08 

0.00 

0.03 

AW1B 

8 

0.00 

1.00 

0.00 

0.02 

0.00 

0.07 

0.00 

SW1T 

9 

0.04 

0.00 

0.95 

0.00 

0.00 

0.00 

0.02 

AW1T 

11 

0.00 

0.04 

0.00 

0.93 

0.01 

0.16 

0.01 

SW2B 

12 

0.07 

0.00 

0.04 

0.00 

0.99 

0.00 

0.34 

AMLGL 

13 

0.00 

0.04 

0.00 

0.31 

0.01 

0.97 

0.01 

SMLGL 

14 

0.03 

0.00 

0.02 

0.00 

0.84 

0.03 

0.61 


*: Nastran mode number 
**: GVT mode number 


Fig. 25 MAC matrix from previous history. 
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Fig. 26 Orthogonality matrix from previous history. 
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Fig. 27 Cross-orthogonality matrix from previous history. 


Table 6 also shows the frequency results for the third model-tuning 
procedure. The MAC matrix after the third model-tuning procedure is 
displayed in Fig. 28. MAC terms for mode number 14 improved for 
both weight configurations. Mass orthogonality matrix results are 
shown in Fig. 29. Small improvements were observed for the EFEW 
configuration with the off-diagonal term 12-13 going from -0. 150 to 
-0.142 and the off-diagonal term 12-14 reducing from 0.107 to 
0.091. Similarly, the off-diagonal term for the FFFW configuration 
went from 0.118 to 0.105. Cross-orthogonality matrix results are 
shown in Fig. 30. Both objective function terms improved, with the 
EFEW mode interchanged off-diagonal term 12-12 improving the 
most from an initial value of 0.238 to 0. 175, whereas the FFFW mode 
interchanged off-diagonal term 14-13 going slightly bad from 0.286 
to 0.288. 

D. Fourth Model-Tuning Procedure 

Based on frequency results after the third model-tuning procedure 
in Table 6, the frequency result for mode number 28 of the EFEW 


configuration stood out due to the fact that frequency error was 
allowed to be at the maximum value of 20%. For this model-tuning 
procedure, the goal was to improve this particular secondary mode 
frequency, and as a result, the frequency of the EFEW mode number 
28 was selected as the objective function for this run. Similar to the 
last procedure, identical design variables were used with constraint 
functions consisting of all other cross-orthogonality matrix terms 
(42 + 42), mass orthogonality matrix terms (21 +21), frequency 
terms (14 + 12), total weight (1 + 1), x c.g. locations (1 + 1), and 
GVT sensor cable weight (1_), for a total of 157 constraint functions. 
The DOT algorithm was used again for this optimization procedure. 

Frequency results for the fourth model-tuning procedure are 
presented in Table 6. Notice that the frequency for the EFEW mode 
number 28, which was the objective function for this run, improved 
from —20% from the previous procedure to — 1 6.7% . MAC results for 
this procedure are shown in Fig. 31_; a slight improvement was 
observed for the FFFW weight configuration, where mode number 
14 went from 0.69 in the last procedure to 0.70 for this run. Figure 32 
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Fig. 28 MAC matrix after the third model-tuning procedure (DOT-03). 
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Fig. 29 Orthogonality matrix after the third model-tuning procedure (DOT-03). 
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Fig. 30 Cross-orthogonality matrix after the third model-tuning procedure (DOT-03). 
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Fig. 31 MAC matrix after the fourth model-tuning procedure (DOT-04). 
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Fig. 32 Orthogonality matrix after the fourth model-tuning procedure (DOT-04). 
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Fig. 33 Cross-orthogonality matrix after the fourth model- tuning procedure (DOT-04). 
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shows the results of the mass orthogonality matrix; only a slight 
improvement was seen for the off-diagonal term 12-13 for the EFEW 
configuration. The off-diagonal term 12-14 for the EFEW and FFFW 
increased from 0.091 to 0. 1 10 and from 0. 105 to 0. 120, respectively. 
Cross-orthogonality matrix results are displayed in Fig. 33. Not much 
change was observed between the last run and this run, with the 
exception of the off-diagonal term 12-14 for the FFFW weight 
configuration, which went from a magnitude of 0.151 to 0.140. 

The design variable changes of the model-tuning procedures are 
summarized in Table 7, which actually displays the design variable 
percentage changes between each model-tuning procedure. Notice 
that the percentage change for many of the main landing gear design 


variables after the fourth model-tuning procedure (DOT-04) 
approached the ±20% side constraint limit. In addition, the smallest 
MAC value, largest frequency error, largest off-diagonal terms of 
orthogonality, and cross-orthogonality matrices were all related to the 
main landing gear. Therefore, it can be concluded that the FE model 
pertaining to the main landing gear was not accurate enough to allow 
significant change or improvement. It was later found that there was 
an idealization error with the main landing gear. This reasoning for 
this idealization error is presented in Fig. 34. In Fig. 34, the main 
landing-gear was originally designed as tapered and curved beams. 
When creating the FE model, the tapered and curved regions of the 
main landing gear were modeled as uniform and straight beams, as 


Table 6 Flexible modes from previous history, and after third and fourth model-tuning procedures 








Nastran results 






GVT data 


POH 



DOT-03 



DOT-04 



Mode 

Mode shape 

Mode 

Frequency 

Error, % 

Mode 

Frequency 

Error 

Mode 

Frequency 

Error 

Target error 

7 

SW1B 

7 

1.086 

1.8 

EFEW configuration 
7 1.090 

2.2 
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1.101 

3.1 

5(3) 
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AW1B 
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0.0 
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1.549 

0.4 
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1.5 

5(3) 
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6.0 

10 

Total weight 


368.1 
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0.3 

-32.19 

— 

z C.g. 
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Table 7 Design variable changes (percent) during model- tuning procedures 



Design variable 

Property 

Base 3 

DOT1 

DOT2 

POH b 

DOT3 

DOT4 

1 

CONM2 c 100202 and 100204 

Mass 

0.0 

0.8 

0.8 

13.4 

17.6 

7.9 

2 

PBAR d 310 

h 

0.0 

0.0 

0.0 

0.0 

0.2 

1.1 

3 

PBAR313 

h 

0.0 

-1.4 

-1.4 

-8.5 

8.9 

11.7 

4 

PBAR314 

h 

0.0 

-0.6 

-0.6 

-3.3 

-10.4 

-20.0 

5 

PBAR 308 

h 

0.0 

-0.5 

-0.5 

-2.0 

-4.7 

-12.9 

6 

PBAR310 

h 

0.0 

0.3 

0.3 

-0.7 

-4.5 

-4.6 

7 

PBAR 313 

h 

0.0 

0.4 

0.4 

0.3 

1.9 

18.0 

8 

PBAR 310 

J 

0.0 

0.0 

0.0 

0.0 

-0.1 

-0.1 

9 

PBAR 313 

J 

0.0 

0.1 

0.1 

0.1 

2.6 

9.2 

10 

PBAR 314 

J 

0.0 

0.4 

0.4 

0.4 

6.4 

20.0 

11 

PBUSH 6 315 and 316 

k 2 

0.0 

0.0 

0.0 

-0.5 

-0.1 

-1.6 

12 

PBUSH 315 and 316 

k 3 

0.0 

0.0 

0.0 

-0.7 

-6.7 

-4.2 

13 

PBUSH 315 and 316 

k 4 

0.0 

0.0 

0.0 

-0.5 

-2.3 

-4.9 

14 

MAT8 f 100 

E 

0.0 

-2.0 

-2.0 

0.7 

1.5 

4.8 

15 

MAT8 100 

G 

0.0 

6.4 

6.4 

9.0 

13.9 

17.9 

16 

PBAR 311 

A 

0.0 

0.0 

0.0 

0.4 

1.0 

1.1 

17 

PBAR 312 

A 

0.0 

0.0 

0.0 

2.3 

9.0 

15.3 

18 

PBAR 313 

A 

0.0 

0.0 

0.0 

1.3 

9.9 

17.1 

19 

PBAR 312 

h 

0.0 

0.0 

0.0 

-0.1 

0.5 

1.0 

20 

PBAR 309 

h 

0.0 

0.0 

0.0 

-1.6 

-6.3 

-17.2 

21 

CONM2 930618 

Mass 

0.0 

0.0 

0.0 

0.0 

-1.9 

-2.2 

22 

CONM2 9120105 

Mass 

0.0 

0.0 

0.0 

1.2 

-20.0 

-20.0 

23 

CONM2 930668 

Mass 

0.0 

0.0 

0.0 

-0.1 

2.0 

2.5 

24 

CONM2 9220105 

Mass 

0.0 

0.0 

0.0 

-2.2 

10.9 

8.1 


a Baseline configuration (tuned by LMSW). 
b Obtained from previous optimization history. 
c Coneentrated mass element (MSC Nastran terminology). 
d Simple beam property (MSC Nastran terminology). 

Generalized spring and damper property (MSC Nastran terminology). 
f Shell element orthotropic material property (MSC Nastran terminology). 


opposed to tapered and curved beams, respectively. Improving 
properties related to the main landing gear was not an easy task 
because only a single uniform sectional property was used for the 
whole tapered region. This was the main reason why some design 
variables in Table 7 stood out due to their large changes, but because a 
±20% side constraint was imposed and none of the design variables 
exceeded this side constraint margin, the resulting changes for the 
physical properties were deemed acceptable. 

Nondimensional flutter speeds and frequencies after the fourth 
model tuning are summarized in Table 3. Corresponding speed 
versus damping V - g and speed versus frequency V - f curves of 
this model with the EFEW configuration at a Mach number of 0.160 
from the matched flutter analyses are given in Fig. 35. Improved 
frequencies and mode shapes in this study increased average flutter 
speeds by 7.6% compared to the baseline model. In Fig. 16, notice 
that LMSW’s model tuning resulted in large increases in flutter speed 
for the baseline model when compared to their nonvalidated final 
design model. From this model-tuning procedure done by NASA, 
further inroads were made by increasing the flutter speed for the 
symmetric body freedom flutter and symmetric wing bending torsion 
flutter, as shown in Fig. _L6. 

VI. Conclusions 

For the X-56A flight-test support, NASA AFRC has undertaken the 
task of improving the GVT test- validated stmctural dynamic finite- 
element model. Using an in-house O 3 tool developed by NASA 
AFRC’s structural dynamics group, along with the Nastran program, 
structural model tuning was conducted so as to improve the correlation 
between numerical and experimental modal data to reduce model 
uncertainties by meeting specified NASA and military standards. 

Stmctural dynamic model tuning of the aircraft was performed in 
stages, starting with tuning of the frequency, followed by the cross- 
orthogonality (mode shape) matrix. A total of four stmctural dynamic 
model tuning mns were conducted focusing on the improvement in 
frequency errors and the off-diagonal terms of the cross-orthogonality 
matrices for both the EFEW and the FFFW configurations of the X-56A 


aircraft in a single optimization mn. Other properties such as weight of 
GVT sensor cables, total weight, c.g. location, and off-diagonal terms of 
the orthogonality matrix were primarily used as constraints because the 
current configurations met the majority of the defined criteria. MAC 
constraints were not used as objective or constraint functions. However, 
the MAC values were used in a mode-tracking program to overcome 
mode interchange problems during model-tuning procedure. 

The first and second model-tuning procedures focused on the 
improvement of frequencies and off-diagonal terms of the cross- 
orthogonality matrices to meet the NASA standard by using main 
landing gear, Young’s modulus, and Shear modulus design variables. 
All of the off-diagonal terms violate military and NASA standards, 
which were subsequently improved after the first two model-tuning 
procedures. An average of 3.8% of improvement was observed in the 
cross-orthogonality matrices. Slight violation of target frequency 
values was also observed after the first two model-tuning procedures. 

Before the third and fourth model-tuning procedures, previous op- 
timization histories were reviewed, and a much better design configu- 
ration was found by change. By selecting an improved optimization 
configuration based on POH as the new starting point, the optimi- 
zation results for the third and fourth model-tuning procedures were 
improved significantly. It should be noted that the largest correlation 
errors were mostly associated with the main landing gear due to the 
existence of the idealization error related to the main landing gear. 

Finally, it should be emphasized that the performance indices 
based on each individual element of frequency errors and the off- 
diagonal terms of the orthogonality and cross-orthogonality matrices 
introduced in this study were easier to use than the norm average- 
based model-tuning method [2,17]. In the case of the norm-based 
approach, the total number of performance indices was a lot smaller 
than the current method; however, it was not easy to control each 
individual off-diagonal term when compared to the current approach. 
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